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The lifetimes for the double beta decays of ^®Ge, *^Se and ^"^^Xe are calculated us- 
ing very large shell model spaces. The two neutrino matrix elements obtained are in 
good agreement with the present experimental data. For < >< 1 eV we pre- 

dict the following upper bounds to the half-lives for the neutrinoless mode: T^'^^!^' (Ge) > 

1.8510^^yr., Tf^l\Se) > 2.3610^^?yr. and Tl'j^\Xe) > 1.2110^^2/r-. These results are 
the first from a new generation of Shell Model calculations reaching 0(10*) dimensions. 
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The importance of the nuclear double beta decay {(3(3) is well established. The two neutrino mode (/3/3(2^) ) that has 
been proven to be very sensitive to the nuclear correlations, provides a severe test of the nuclear wave functions. The 
neutrinoless mode (/3/3(oiy)) is one of the best probes for the physics beyond the Standard Model. It is particularly 
suitable to explore the intrinsic properties of the neutrino like its mass and the existence of right-handed weak currents 




Reasonable values of these lepton number violating parameters can be extracted from the experiment provided that 
accurate nuclear matrix elements are used. Previous to this work, large Shell Model (SM) calculations of the P(3 decay 
were only possible in '^^Ca [0J^,^. Then, several approximations had to be used in order to study heavier nuclei. 

A weak coupling approximation was used by Haxton, Stephenson and Strottman in order to calculate the decays 
of ''^Ge, ®^Se and i^s.iaorpj^ They also invoked the closure approximation to circumvent the calculation of the 1+ 
states in the intermediate nuclei. Therefore, it is very difficult to estimate the uncertainty in their results because 
there is not an unique method to choose the energy denominator of the two neutrino matrix element. Using an 
statistical method to determine the energy denominators and, setting qa = 1, they gave reasonable values for the 
lifetimes in all the cases, except for the Tellurium isotopes. 

The quasiparticle random-phase approximation (QRPA) has also been used to study the decay of medium and heavy 
nuclei . It was found ||^,^ that the matrix elements of the Gamow- Teller and double Gamow- Teller operators were 
very sensitive to the particle-particle interaction. The strength of this interaction (gpp) is treated as a parameter and 
fitted to the available /3± data. The introduction of gpp as a phenomenological factor makes it possible to reproduce 
the observed half-lives. However, the predictivity of the QRPA approach is severely limited because of the large 
variation of the relevant matrix element in the physical window for gpp. The neutrinoless mode is less sensitive 

to the correlations, but there still exists a factor of three among several predictions for Mq^'^ ||1c|-[T^ 

Very recently the two neutrino mode has been studied in the Shell Model Monte Carlo method (SMMC) For 
the decay ^^Ca Ti the matrix element is 0.15±0.07 MeV~^ Using the same interaction, the exact diagonalization 
gives 0.08 MeV-^ For "^^Ge a matrix element 0.13 ± 0.05 MeV"^ is found. 

The aim of this paper is the study of the /3/3 decay of ^^Ge, ®^Se and ^^^Xe in the SM framework. This study presents 
a new generation of SM calculations that have recently become feasible due to the availability of a very performing 
SM code [l^ that makes it possible to carry the calculations in an ordinary workstation. Similar importance has 
the development of a minimal phenomenological treatment of the realistic interactions in order to cure their bad 
saturation properties jisj. We believe that these calculations represent a real improvement over the old SM and 
QRPA calculations 

For ''^Ge and *^Se the Schrodinger equation is solved in a valence space that consists of the following orbitals: P3/2, 
/5/2J Pi/2 and (79/2. ^^^Xe is studied in a valence space that is made of the ^5/2, Si/2, <?7/2,'^3/2 and ft.11/2 shells. 
In table | we show the number of active particles as well as the dimensions reached in each case. In the first space 
t measures the number of particles that are allowed to jump from the pf subshells to the gg/2 orbital. In the case 
of A=136 it represents the number of particles that can be excited from the lower (77/2 and subshells to the 
remaining three orbitals. In both cases "full" indicates that no limitations are imposed. Later on we shall analyse the 
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convergence of the relevant matrix elements from t=0 to the full calculation, except in A=76 where we are limited to 
t=4(Ge and Se) and t=5(As) because the full calculation is out - but not far - of reach. 

We use as starting point a G-matrix from Kuo to calculate the wave functions of the A=76,82 nuclei . To fix 
the interaction the monopole parameters are fitted to the energy levels of the Ni isotopes and the N=50 isotones [ pT| . 
The main monopole changes amount to weaken the interaction among the pf orbits and the gg/2 shell. For ^^^Xe, 
the interaction is the G-matrix obtained from the the Bonn potential with monopole modifications in order to 
reproduce the spectroscopy of the N=82 isotones The secular problem is solved using the Lanczos algorithm. 

The eigenstates are used for calculating the matrix elements of the weak interaction inducing the /3/3 decays. We 
describe the weak processes in the nucleus by an effective hamiltonian proposed by Doi et al. ^ , that consists of V 
and A currents and that is compatible with SU{2)l <Xi SU{2)fj (g) C/(l) grand unification models. 
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In addition to the right-left and right-right coupling constants, the neutrino mass (implicit in the leptonic currents) 
also breaks the maximal parity violation of the standard theory. This hamiltonian led us to an scenario where there 
is no maximal parity violation and where the leptonic number is not preserved. The neutrinoless double beta decay is 
mediated by (virtual) massive Majorana neutrinos. Moreover the description of the two neutrino mode is essentially 
the same as that of the Standard Model. This is so because the relevant contributions come from the standard left-left 
coupling. With regard to the hadronic currents, all the terms up to order v/c are included because, as first noticed 
by Tomoda et al. [ p^ , they contribute significantly to the zero neutrino mode. 
(/3/9)2i/ results: The half- life can be approximated as 
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where G is a Integral kinematical factor and M^!^ is the usual energy weighted double Gamow- Teller matrix element. 
We use the bare Gamow- Teller operator in the definition of the matrix element and use = 1 in the integral kinematic 
factor. The exact definition of this matrix element as well as the algorithm that we use to calculate it can be found 
in reference Q . By means of this algorithm we can get a reliable approximation to the Gamow- Teller strength in the 
intermediate nucleus. However, due to the fact that our model spaces do not include all the spin partner orbits we 



can not exhaust the 3(N-Z) sum rule. In fact, we find: 5 



(Ge) 



17.14, S. 



(Se) 



21.66 and S". 



(Xe 



52.30 compared 



to the 3(N-Z) values 36, 42 and 84 respectively. However, as the remaining strength will appear at high energy its 
influence in the description of the two neutrino mode should be small. 

We summarize our results in table |l[ This is the only mode whose existence has been clearly confirmed Q and 
the lifetimes have been measured for several emitters. There are two sets of matrix elements. In the first one the 
theoretical energies of the 1+ states are used whilst, in the second the spectrum of these states is globally shifted in 
order to place the first 1+ state at its experimental energy. As it can be seen in table || the matrix elements increase a 
20% when the experimental energies are used. The agreement with the experimental data is, in both cases, reasonably 
good. Further refinements of the interaction are possible but, we do not expect major modifications of the matrix 
elements. 

It is worth comment that our closure matrix elements are very different from those of reference ||l[, 0.68 compared 
to 2.56 (''^Ge) and 0.74 compared to 1.876 (^^Se). These discrepancies also led to quite different effective 1+ centroids. 

In the A=76,82 region the matrix element increases very slowly as the valence space is enlarged. Consequently, the 
half life decreases and there is a factor two between the t=2 and the final predictions. This behavior is different to 
that found in the decay of ^^^Xe, where the matrix element is nearly constant. 

It is possible to compare the SM result, 0.14 with the SMMC extrapolation, 0.13±0.05. Both numbers are very close. 
Nevertheless, the comparison of other relevant quantities {{M^'^j,)c, S+, E) is not as good. Any further discussion 
must be delayed because, although the valence spaces are identical, the effective interactions are different. 

(/3/9)oi/ results: In the closure approximation the half- life of the 0+— > 0+ decay can be written as 
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where < >, < A > and < 77 > are the effective lepton violating parameters, V'i(2) are the CP phases and the Cxy 
coefScients are hnear combinations of the nine matrix elements and nine integral factors. A clear and comprehensive 
definition of them all can be found in reference 

The nine matrix elements shown in table [11 are calculated in the light neutrino approximation (see Doi |^) and 
using qaI gv = 1-25. Since this mode has not been observed yet, we take a half-life of lO^^yr., which is very close to 
the expected experimental limits . This makes it possible to obtain the upper bounds to the three lepton violating 
parameters < >, < A > and < rj >. The results are compiled in table IV. For < rriiy >< leV we find the following 
lower bounds: T^^/2^(G'e) > LSSlO^^yr. Tl°^\Se) > 2.36102V- and T^^{Xe) > 1.21 lO^^yr. Although the result 
for ^^Ge has been obtained in a t=4 calculation, we expect that the matrix element and therefore the half-life are close 
to convergence, as it is the case in ^'^Se. Being the {I3(3)2v mode a background for the detection of the neutrinoless 
mode the ratio T^^yj very important. Notice that, among the three studied nuclei, this ratio is minimum for 

^"^Se and, therefore, it can be considered a good candidate for the detection of the zero neutrino mode. 

There are discrepancies among different calculations of the I3(3{pu) rnode. Table ^ compares the Mq!^ values 
predicted by several authors. We see that our matrix elements are smaller than the others. The only exception are 
the new results by Pantis et al. , that include to some extent the neutron-proton pairing interaction. Based on the 
reasonable agreement between theory and experiment for the spectroscopy of the regions around these (3(3 emitters 
and for the I3[3{2v) mode, we are also quite confident on our predictions for the zero neutrino mode. 

In summary, we have presented the first large scale SM calculations of the (3(3 decays of ^^Ge, ^^Se and ^■^^Xe. The 
agreement with the present experimental data, available for the I3[3{2v) mode, is reasonable and encouraging. For the 
f3P(^Q,^) mode we have obtained upper bounds to the neutrino mass and to the coupling constants assuming a lifetime 
of lO^^yr., and the lifetimes predictions for < >< leV. 

Acknowledgments. We thank M. Hjorth-Jcnscn and T.T. Kuo for making their G-matrix available to us. We 
thank D. Dean, P. Vogel and A. Zukcr for their help. This work has been partly supported by the IN2P3 (France)- 
CICYT(Spain) agreements and by a grant of DGICYT (Spain) PB93-263. 



W. C. Haxton and G. J. Stephenson, Jr., Progress in Particle and Nuclear Physics 12, 409 (1984). 
M. Doi, T. Kotani and E. Takasugi, Progress of Theoretical Physics Supplement 83, 1 (1985). 
F. Boehm and P. Vogel, Phystcs of Massive Neutrinos, Cambridge: Cambridge Univ. Press. (1992). 
M. Moe and P. Vogel, Ann. Rev. Nucl. Part. Sci. 44, 247 (1994). 

E. Caurier, A. Poves and A.P. Zuker, Phys. Lett. B 252, 13 (1990). See the erratum in E. Caurier, A. P. Zuker, A. Poves 

and G. Martmez-Pinedo, Phys. Rev. C50, 225 (1994). 

J. Retamosa, E. Caurier and F. Nowacki, Phys. Rev. C 51, 371 (1995). 

P. Vogel and M. R. Zimbauer, Phys. Rev. Lett. 57, 3148 (1986). 

O. Civitarese, A. Faessler, T. Tomoda, Phys. Lett. B 194, 11 (1987). 

K. Muto and H.V. Klapdor, Phys. Lett. B 201, 420 (1988); K. Muto, E. Bender and H. V. Klapdor, Z. Phys. A - Atomic 
Nuclei 334, 177 (1989). 

T. Tomoda, A. Faessler, Phys. Lett. B 199, 475 (1987). 

K. Muto, E. Bender and H. V. Klapdor, Z. Phys. A - Atomic Nuclei 334, 187 (1989). 

A. Staudt, T. T. S. Kuo and H. V. Klapdor , Phys. Lett. B 242, 17 (1990). Atomic Nuclei 334, 187 (1989). 
P. B. Radha et al, Preprint |nucl-th9510052 
E. 
A. 
T. 
F. 
M 
T. 



Caurier, SM code ANTOINE, ( 1989). 

Zuker and M. Dufour, Preprint ^lucl-th9505012| , Preprint |nucl-th950501C 
T. Kuo, private communication 



Nowacki, Phd thesis, ULP Strasbourg (1995), unpublished 
Hjorth- Jensen, T. T. S. Kuo and E. Osnes; Phys. Rep. 261, 125 (1995). 
Tomoda, A. Faessler, K. W. Schmid and F. Grummer, Nucl. Phys. A452, 591 (1986). 
NEMO 3 Proposal, LAL 94-29 (1994). 

H. V. Klapdor-Kleingrothaus, Proc. Inter. Workshop on Double Beta Decay and Related Topics (Trento 1995). Eds. H. V. 

Klapdor-Kleingrothaus and S. Stoica (World Scientific Singapore). 

S. R. Elliot et al. Phys. Rev. C46, 1535 (1992). 

K. Lou et al., Proc. of Osaka Conference, June 1995. 

G. Pantis, F. Simkovic, J. D. Vergados and Amand Faessler; to be published. 



4 



TABLE I. Dimensions and active nucleons for each nucleus. See text for the meaning of t. 









Dimensions 


Truncation 




4 


16 


8176629 


t = A 




5 


15 


63788368 


t = 5 




6 


14 


54625321 


t = 4 


82 Op 


u 






full 




7 


19 


8353667 


full 




8 


18 


70757366 


full 




4 


32 


1504 


full 




5 


31 


212338 


full 


i36Ba 


6 


30 


13139846 


full 



TABLE II. Two neutrino matrix elements and half lives. Mq,^^ in MeV ^ and T1/2 in years. 



Decay t (M^^p rff^ {Tlj^)e.p (M^t^)c S. 

AE{exp) AE{th) AE{exp) AE{th) 










0, 


.000 


0, 


.000 














0, 


.000 


16, 


.73 


0.000 


^^Ge-. 




2 


0, 


.112 


0, 


.088 




5, 


.6781021 


9, 


.1971021 




0, 


.465 


17, 


.13 


0.146 






4 


0, 


.180 


0, 


.140 




2 


.1981021 


3, 


.6341021 




0, 


.676 


17, 


.14 


0.258 
















0.22 










1.80 1021 



















0, 


.000 


0, 


.000 














0, 


.000 


21, 


.66 


0.000 




82 


2 


0, 


.128 


0, 


.102 




1 


.312 1020 


2, 


.065 1020 




0, 


.483 


21, 


.61 


0.121 






4 


0, 


.198 


0, 


.155 




5, 


.482 IQi^ 


8, 


.946 IO19 




0, 


.745 


21, 


.56 


0.209 








0, 


.208 


0, 


.164 


0.14 


4, 


.968 IQi^ 


7, 


.991 IO19 


1.08 1020 


0, 


.799 


21, 


.55 


0.226 









0, 


.026 


0, 


.028 




2, 


.4871021 


2, 


.455 1021 




0, 


.106 


52, 


.75 


0.004 


i36Xe - 




2 


0, 


.036 


0, 


.039 




1 


.485 1021 


1, 


.265 1021 




0, 


.178 


52, 


,37 


0.007 






4 


0, 


.032 


0, 


.035 




1 


.8791021 


1, 


.571 1021 




0, 


.146 


52, 


.30 


0.008 






full 


0, 


.031 


0, 


.034 


< 0.06 


2, 


.003 1021 


1, 


.665 1021 


> 5.60 1020 


0, 


.143 


52, 


.30 


0.008 
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TABLE III. Nuclear matrix elements for the (/?/?) (o,.) mode. Xqj. = 2 - x'qj., Xp = 2Xf - x'p. 



M.E. 






= 


t 


= 2 


t ■ 


= 4 


Full 






0, 


.721 


1, 


.294 


1, 


.568 










0, 


.505 


1, 


,340 


1, 


,846 


1 

1 , 


.y * u 






0, 


.484 


0, 


.630 


0, 


.649 


u. 


.001 


XF 


TP. y' 1 


-0, 


.068 


-0, 


.098 


-0, 


.106 








S9 n 


-0, 


.101 


-0, 


.107 


-0, 


.107 


— u 


1 HQ 
.iUo 






-0, 


.172 


-0, 


.156 


-0, 


.157 


— u 


.ioo 


Xgt 


"^^Ge 


1, 


.074 


1, 


.107 


1, 


,115 








^^Se 


1, 


.069 


1, 


.114 


1, 


,119 


i, 


.120 




^^^Xe 


1, 


.103 


1, 


.106 


1, 


.099 


i, 


.097 


Xf' 


"^^Ge 


-0, 


.060 


-0, 


.102 


-0, 


.109 










-0, 


.104 


-0, 


.112 


-0, 


.112 


— 


.112 






-0, 


.184 


-0, 


.166 


-0, 


,167 


— 


.167 


XT 


^^Ge 


0, 


.186 


0. 


.043 


0, 


.017 










0, 


.156 


0, 


.049 


0, 


.031 


0, 


.028 






0- 


.039 


-0. 


.006 


-0, 


.031 


-0 


.031 


y p 


76 Ge 


-1 


.435 


-0, 


.832 


-0, 


.544 










1, 


.710 


0, 


,852 


0, 


,574 


0, 


.494 






0, 


.898 


0, 


.411 


0, 


,280 


0, 


.256 


Xr 


76Ge 


0, 


.761 


0, 


.707 


0, 


.684 










0, 


.873 


0, 


.706 


0, 


.683 


0, 


.680 




i36Xe 


0, 


.780 


0, 


.872 


0, 


.942 


0, 


.955 



TABLE IV. Shell Model predictions for the upper bounds to the three lepton violating parameters (LPV): the effective mass 
< rriv > and the coupling constants < A > and < r; > for T1/2 > W^^yr. 









t = 


t = 2 


t = 4 


Full 


76Ge 


< 


> (eV) 


2.98 


1.60 


1.32 






< 


rj > 10« 


2.78 


1.83 


1.65 






< 


A > 10^ 


4.31 


2.65 


2.24 






< 


m„ > (eV) 


1.93 


0.72 


0.52 


0.49 




< 


rj > 10« 


3.03 


1.22 


0.86 


0.79 




< 


A > 10^ 


1.93 


0.82 


0.61 


0.57 




< 


m„ > {eV) 


1.47 


1.14 


1.10 


1.10 




< 


r] > 10« 


3.10 


1.77 


1.53 


1.49 




< 


A > 10^ 


1.88 


1.53 


1.50 


1.49 
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TABLE V. SM versus QRPA M^^' matrix elements. 





This Work 


Haxton 


1 Tomoda 


Muto |Tl| 


Staut H 


Pantis Hi 




1.568 


4.180 


3.355 


3.014 


10.910 


1.846 




1.970 


3.450 


3.055 


2.847 




1.153 




0.651 




1.120 






1.346 
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